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^ — ' . Abstract. The narrow K-z relation of powerful radio galaxies in the Hubble K diagram is often attributed to the 

' stellar populations of massive elliptical galaxies. Because it extends over a large range of redshifts {0 < z <4), 

04 I it is difficult to estimate masses at high redshifts by taking into account galaxy evolution. In the present paper, 

we propose to estimate the stellar masses of galaxies using the galaxy evolution model PEGASE. We use star 
formation scenarios that successfully fit faint galaxy counts as well as z—0 galaxy templates. These scenarios also 
predict spectra at higher z, used to estimate valid photometric redshifts. The baryonic mass of the initial gas 
cloud Mbar is then derived. The K-z relation is remarkably reproduced by our evolutionary scenario for elliptical 
galaxies of baryonic mass Mbar, max —10^^ Mq, at all z up to z= 4. Mbar, max is also the maximum mass limit of all 
types of galaxies. Using another initial mass function (IMF), even a top-heavy one, does not alter our conclusions. 
The high value of Mbar, max observed at z > 4 implies that massive clouds were already formed at early epochs. 
fT^ ] We also find that the Mbar, max limit is similar to the critical mass Afcrit of a self-gravitating cloud regulated by 

. cooling (Rees & Ostriker, 1977; Silk, 1977). Moreover, the critical size rcrit — 75 Kpc is remarkably close to the 

typical diameter of Lya haloes surrounding distant radio galaxies. This confirms the validity of the method of 
baryonic mass determination based on the K-hand luminosity. A puzzling question that remains to be answered 
is the short time-scale of mass-accumulation required to form such massive galaxies at z = 4. We discuss the 
O ' dispersion of the K-z relation in terms of uncertainties on the mass limit. The link between the presence of the 

, active nucleus and a large stellar mass is also discussed. 
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1. Introduction because galaxies are more gas-rich and formed stars more 

Ci actively in the past. 
■ The K Hubble diagram is an efhcient tool to study the 

evolution of galaxies at high redshifts. Because evolution- The main feature of the galaxy distribution in the 

ary processes in the expanding universe are not known, Hubble A'-band diagram is the well-known K-z relation, 

models are required. The most valuable models propose the sharp bright limit traced by nearby (z ~ 0) massive 

scenarios of star formation that aim to fit observations of galaxies and by hosts of distant powerful radio galaxies. 

galaxies at all redshifts; the first constraint is to reproduce The physical meaning of the relation, up to z= 4 or more, 

the templates of nearby (z ~0) galaxies. is still puzzling. Most of the interpretations of the K-z 

The K-hand (centred at 2.2 ^m) is preferred for evo- relation explain its dispersion either by testing various 

lutionary analysis because galaxies are seen over a large sets of cosmological parameters (Hq, qo, Aq) (Longair & 

redshift range (0< z < 4) in their near-IR to optical rest Lilly, 1984, and more recently Inskip et al. 2002) or by 

frames. Moreover, the evolutionary effects are reduced, examining the relation of the galaxy K magnitude with 

The Hubble K-z diagram is known to be an excellent tool its radio power. The radio-powerful SCR galaxies are on 

to measure stellar masses of galaxies, up to high redshifts. average brighter {AK ~ - 0.6 mag) than the less radio- 

At z ~0, the K luminosities are dominated by the bulk of powerful 6C (Bales & Rawhngs, 1996; Bales et al, 1997) 

old low mass red stars. At higher z, the K emission is due or 7C (Willott et al. 2003) radio galaxies. Moreover the 

to the redshifted emission of blue luminous young stars radio-powerful Molonglo sample (McCarthy, 1999) shows 

a AK ~ — 1 mag with the 6C catalogue at all redshifts. 
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Fig. 1. Radio galaxy hosts (empty squares) and deep field galaxies (dots and crosses) plotted in the Hubble K 
diagram. The catalogues of radio galaxies are respectively SCR (blue squares) and 6CE (green squares). The sample 
also includes highly distant radiogalaxies observed with the NIRC and some 6C* radio galaxies (red squares). The 
field galaxy samples are the HDF-N (black crosses, Williams et al., 1996) and Hawaii (black dots, Cowie et al., 1999) 
surveys. 



lations of radio galaxy hosts are identified with massive 
elliptical galaxies, even at high redshifts (van Breugel et 
al., 1998; Lacy et al., 2000; Pentericci et al., 2001). Radio 
galaxy luminosities are generally high, typically L ~ 3 
to 5 L, for both radio-loud and radio-quiet galaxy hosts 
(Papovich et al., 2001). Radio galaxy hosts are located 
in place of elliptical galaxies in the fundamental plane 
(Kukula et al., 2002). At high z, mass estimates of ra- 
dio galaxy hosts are rare. From an incomplete rotation 
curve, Dey et al., (1996) derived a mass of 10^^ for 
3C 265 (z=0.81). Using HST/NICMOS observations, Zirm 
et al., (2003) confirmed that radio galaxy hosts are mas- 
sive, even at high z. Stellar masses estimated by fitting 
stellar energy distributions (SEDs) to the NICMOS data 
are between 3 and 8 10" Mg. The CO (3-2) fine of the 
very massive galaxy SMM J02399-0136 (z=2.8) observed 
with the IRAM observatory by Genzel et al., (2003) shows 
that it is a rapidly rotating disk with a total dynamical 
mass of ~ 3 10^^ sin^ i M© {i is the galaxy inclination). 
Stellar masses of three galaxies, selected in the Hubble 
Deep Field South at Kg < 22 on the basis of their un- 
usually red near-IR color (J-K > 3), were estimated to be 
about 10" Mq from their SEDs (Saracco et al., 2003). 

Section 2 presents the observational samples of radio 
and field galaxies in the Hubble K diagram. Following 
a preliminary interpretation of the K-z relation (Rocca- 
Volmerange, 2001), we explore in section 3 the space of 
cosmology and galaxy evolution parameters with the code 
PEGASE to confirm this interpretation. We successively 
analyse the parameter set: distance modulus, star forma- 



tion scenarios that define the evolutionary time-scales of 
the various components, the IMF and finally the baryonic 
mass of the initial gas reservoir. For the sake of clarity, evo- 
lution scenarios are limited to spiral and elliptical types: 
considering scenarios between these two types would not 
modify our conclusions. In section 4, the K-z diagram is 
analysed with the predicted K-z sequences of spiral and 
elliptical galaxies formed from clouds of various baryonic 
masses Mbar=10^ to 10^^ M©. We find that only the se- 
quence of the elliptical galaxies formed from a 10^^ Mq 
progenitor cloud uniformly fits the K-z relation from z=0 
to z=4. The sensitivity of this result to other factors (cos- 
mology models, intensities of emission lines) is considered. 

The striking similarity of the Afbar,max=10^^ Mg mass 
limit to the theoretical estimate of the critical mass of 
fragmentation for self-gravitating clouds (Rees & Ostriker, 
1977) is finally discussed in section 5. Since the K-z rela- 
tion is mainly traced by the most powerful radio galaxies, 
this means that the most massive galaxies correspond to 
the most massive black holes. Masses of galaxies of var- 
ious radio powers are derived in section 6. In particular, 
strong constraints on mass accumulation time-scales are 
given by the high masses of radio galaxy hosts at z ~4. 
The discussion and conclusion make up the last sections. 

2. The observational K-z relation 

The observational sample of distant field galaxies and of 
radio-loud galaxies in the X-band has been gathered by 
De Breuck et al., (2002) (Fig. Radio sources are from 
the 3CR and 6CE catalogues completed by radio galaxies 
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Fig. 2. (a) Plot of the absolute bolometric magnitude Mboi(z, t(z)) vs z for starburst, spiral Sc and elliptical template 
galaxies, (b) similar plot with the apparent bolometric magnitude mboi(z, t(z)). The comparison of the two plots shows 
that the distance modulus is the primary explanation for the K-z slope up to z ~ 1, whatever the spectral type. At 
higher z, stellar evolution has a major effect on light emission for all spectral types. 



observed with NIRC at the Keck I telescope (Laing et al., 
1983, Bales et al., 1997, van Breugel et al., 1998). All the 
magnitudes of radio sources are corrected to a standard 
64 kpc metric aperture, following the method of Bales et 
al., (1997). The brightest {z ~1) SCR sample is complete. 
The NIRC sample is a selection by van Breugel et al., 
(1998) of powerful radio sources at the highest redshifts 
(1.8 < z < 4.4) with strong emission lines. The 6CB and 
NIRC catalogues are not complete. In sections 4 and af- 
ter, the complete 7C catalogues and the 6CE* sample, 
recently published by Willott et al. , (2003) , were added to 
the De Breuck's sample for the interpretation of the K-z 
relation. Over the 8-magnitude range of i^-band luminosi- 
ties of galaxies, the bright envelope of the radio galaxy 
distribution delineates a sharp limit, up to z=4: the so- 
called relation. For comparison, field galaxies of the 
Hawaii survey (Cowie et al., 1999) and the HDF-N survey 
(WiUiams et al., 1996), observed in the K-h&nd at KPNO 
by Dickinson et al., (2003) are also plotted in the Hubble 
diagram. The spectroscopic redshifts were published by 
Fernandez-Soto et al., (1999), Cohen et al., (2000) and 
Songaila et al., (1994). 



3. Modeling galaxy evolution through the K-hanA 

The variation of apparent K magnitudes with redshift is 
predicted by galaxy evolution models, assuming a cos- 
mological model. To clarify the sensitivity of the model 
outputs to the adopted star formation scenarios, we sep- 
arately analyze the effects of other parameters (distance 
modulus, IMF, initial mass and redshift, cosmological pa- 
rameters). Stellar energy distributions (SED) of galaxy 
templates are predicted at all redshifts. The main emission 
lines, due to ionization by massive stars, are computed 
with an HII region model in which the Lyman contin- 
uum photon number and metallicity are the main param- 



eters. Their intensities are added to the SEDs, assuming 
a FWHM of 10 A at z=0. 

For each spectral type, the reddened apparent magni- 
tude of a synthetic galaxy through the filter A at redshift 
z and cosmic time t{z) is given by: 

mx{z,t{z)) =A/A(0,i(0)) 

-f (to - A/)boi(2) + Ax{z) + kx{z) + ex{z) 

where M;)j(0, t(0)) is the intrinsic magnitude at z = and 
at the present cosmic time t(0). The distance modulus 
(to — M)f,o;(z) will be considered below. The k-correction 
k\{z) and e-correction e\(z) respectively depend on the 
cosmology and on the evolution scenario. The extinction 
A\{z) through the filter A is derived from the total amount 
of metal and dust geometry. In our code PEGASB, the 
extinction factor for elliptical galaxies is computed with 
a radiative transfer code for a dust distribution fitted on 
a King's profile and Draine's grain model (see Fioc, 1997 
for details). For spirals, the stars, gas and dust are mixed 
homogeneously. 

We adopt the cosmological model of Friedmann- 
Lemaitre with the Hubble constant Hq in km. s~^. Mpc~^ 
and the matter and cosmological constant den- 
sity parameters. We successively consider the Hq, ^Im, 
values: (65, 0.3, 0.7) for standard , (65, 0.1, 0.) for open 
and (65, 1., 0.) for flat universes. Note that evolutionary 
models are able to successfully reproduce a wide range 
of galaxy observables only in the open and non-null f^A 
universes. 

3.1. The distance modulus 

The distance modulus (m — M)boi(2:) is related to the lu- 
minosity distance Di^{z), and depends only on the cos- 
mological parameters (Weinberg, 1972). Fig. |2] shows the 
absolute (a) to apparent (b) magnitude change in the 
standard cosmology for three spectral types (starburst. 
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elliptical, spiral Sc). At low z, the evolutionary effects of 
early-type galaxies are negligible because the emission is 
dominated by the bulk of low-mass giant stars. The slight 
luminosity increase for spirals is due to recent star forma- 
tion. At higher z, evolution effects become dominant when 
star formation rates are more intense, in particular when 
z ~ Zfor (=10 in our model); the galaxy formation redshift 
Zfor is the redshift when simultaneously infall and initial 
star formation begin. 

3.2. Star formation rates and mass fractions 

In the following, Afbar is the initial mass of the gas reser- 
voir, the baryonic mass of the progenitor cloud. The sce- 
narios of star formation conservatively assume star forma- 
tion rates proportional to the current gas mass. The dura- 
tion of star formation decreases from irregular and spiral 
(> 10 Gyr) to elliptical (< IGyr ) galaxies. The classical 
formalism is only briefly reviewed here. The number of 
stars formed per log mass unit and time unit is 

d?N{m,t) = SFR(t) X IMF(m)dlog(m)dt. 

SFR(t) — Mgas{t)/Tsf is the star formation law of the 
galaxy where Mgas [t) is the current gas mass of the galaxy 
and Tsf is the star formation efflciency. IMF is the initial 
mass function and is assumed not to evolve. All variables 
are normalized to the mass unit of the initial cloud, so that 
the model outputs are multispectral luminosities (L/M)^ 
and mass fractions (stellar: Mgtars, gaseous: Mgas). More 
details are given in the README of the code PEGASE.2 ^ 
(Fioc & Rocca-Volmerange, 1997, 1999c). For simplicity, 
among the eight scenarios of galaxy evolution proposed 
by PEGASE.2, we only select those for elliptical and spi- 
ral Sc galaxies. The extreme difference between their star 
formation efflciencies Tsf implies different star formation 
durations (~ IGyr for ellipticals and > 10 Gyr for spi- 
rals). This is long enough to constrain the accumulation 
time-scales of the stellar and progenitor masses. Another 
reason to adopt the two scenarios is that the formation law 
of our spiral Sc scenario is close to the star formation law 
deduced from the hierarchical model by Baron & White 
(1987) while our law for ellipticals is drastically different 
to hierarchical predictions. Note that the instantaneous 
starburst, often adopted to model stellar populations of 
radio galaxy hosts (Ridgway et al., 2001, Willott et al., 
2003), is unphysical. Moreover, because the time-scale is 
too short, the scenario gives too red colors when they are 
compared to those of massive elliptical galaxies, typical of 
powerful radio galaxy hosts. 

Fig. O presents, respectively for elliptical (a) and spi- 
ral (b) galaxies, the evolution with age and z of the 
galaxy mass Mgai and of its stellar Mgtars and gaseous 
Mgas masses. For the two types, star formation starts at 
infall from the initial gas reservoir. It simultaneously initi- 
ates galaxy formation, so that the initial galaxy mass Mgai 
is null. Mgas measures the gas mass within the galaxy. 

^ Available at ^lLttp://www.iap.fr/pegase| 



Neutral atomic HI and molecular H2 hydrogen are not 
distinguished. The stellar mass Mgtars only corresponds to 
stars still alive while stellar remnant masses are taken into 

account in Mrcmnants- Mgai is the sum of Mstars + Mgas + 

Mromnants- Template elliptical scenarios are characterized 
by a decreasing infall time scale of 300 Myr and the star 
formation efflciency Ts/=0.3 Myr.Mg^. Then supernovae 
produce strong galactic winds starting at 1 Gyr, eliminat- 
ing all gaseous components and instantaneously stopping 
any further star formation. Note that our scenarios do not 
assume further inflow of cool gas from the halo. From Fig. 
in^t), the elliptical scenario is then characterized by peaks 
of Mgai and Mstars as Mgas drops because galactic winds. 
After the age of IGyr, the difference between Mgai and 
Mstars of ellipticals is due to the integrated mass of stel- 
lar remnants Mj-cmnants- In template Sc spirals, the infall 
time scale is 8 Gyr and the efficiency is Tsf= lOMyr.MQ^ 
without galactic winds. The star formation time duration 
is much longer (> 10 Gyr). The baryonic mass M^ar of the 
progenitor is derived by fitting the apparent K luminosi- 
ties Lk of galaxies with the predictions (L/M)k x M^ar 

for IQS < Mfcar/Mo < 1013. 

3.3. Initial mass function 

The IMF parameters are the slope parameter x 
{dN{m) /dlogm a m^) and the lower and upper mass lim- 
its. Three IMFs are considered. First, for Rana & Basu's 
(1992) IMF, the slope parameter above 1.4 Mq is x—- 
1.51 and X—-1.71 above 6.5 M0. Second, for Salpeter's 
(1959) IMF, the slope parameter is x— -1.35. Finally, for 
Gibson & Matteucci's (1997) IMF, x= -0.8, which greatly 
increases the mass of metal ejected. This " top-heavy" IMF 
was proposed by the authors to explain high metal abun- 
dances at the most remote epochs. For all IMFs, we adopt 
a lower mass of 0.1 Mq and an upper mass of 120 M©. 

A truncated IMF, biased towards massive stars, with 
a lower mass of 1 Mq and a slope of -0.8 (the upper mass 
is unchanged) is unable to reproduce the bright iiT-band 
magnitudes of nearby elliptical galaxies: the bulk of low 
mass stars that dominates the stellar emission at z=0 in 
the K-hand is then missing. It will not be considered be- 
low. 

4. The K-z relation 

Figure ^ compares the observed distributions of galaxies 
in the Hubble K diagram with the predictions of elliptical 
(a) and spiral (b) models. The model sequences are com- 
puted with constant values of Mbar in the range 10^ to 
10^3 Mq for the two scenarios. The IMF is from Rana & 
Basu (1992) for all sequences, except for the IO^'^Mq se- 
quence also computed with the "top-heavy" (TH) IMF. 

All observations are covered by the sequences of Mbar 
between 10^ and 10^^ M0; so the initial mass of the reser- 
voir cloud Mbar is the main parameter explaining the evo- 
lution of stellar masses in the Hubble isT-band diagram. 
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Fig. 3. Evolution of mass fractions (Mgtars, Mg^s and Afgai=Afstars+-^^gas+ Afrcmnants) for elliptical (a) and spiral (b) 
galaxies. Ages in Gyr and corresponding redshifts for a standard cosmology are indicated. The formation redshift is 
Zfoi— 10. All masses vary between to 1, after dividing by the initial gas cloud of mass Mbai - Peaks observed at about 
1 Gyr in the elliptical plots result from the stopping of star formation activity by galactic winds following the infall 
phase. This scenario fits the template SED of z ~ elliptical galaxies (see text for details). 



Fig. 0^ also shows the continuous fit of the bright 
K-z relation by the sequence of elliptical galaxies of 
Mbar,max 10^^ M©, up to z > 4. For Spirals, the 10^^ 
sequence (Fig. 01d) is not acceptable at high z while the 
10^'^ Mq one, possible at high z, has no match at low z. 

Confirming that powerful radio galaxies are hosted by 
massive ellipticals (van Breugel et al. 1998; De Breuck et 
al. 2002; Willott et al. 2003), we also estimate their stellar 
Mstars and galaxy Mgai masses at all redshifts: for ellipti- 
cals at z=0 (14 Gyrs) (Fig.|3i) Mstars ^ 55% and A/gai ~ 
65% of the initial mass Afbar of the gaseous progenitor; 
these values only vary within ±10% between z—A and z=Q. 
Thus ~ 35% of the baryonic mass is still in the halo and 
Afromnants — 10%. Gas cxchangcs as late infiows may oc- 
casionally change these fractions. However, depending on 
their efficiency, they would induce fluctuations of the star 



formation history incompatible with the uniformity of the 
K-z relation. 

The scenario of spirals (Fig. ^p) gives, at z={) (14 
Gyrs), A/stars ^ 35%, Afgas ^ 45% and A/gai ~ 85% of 
the initial mass Afbar- These masses rapidly evolve with 
redshift: at z ~2, Afstars ^ 3%, Afgas ^ 27% and A/gai ~ 
31% of the initial mass Afbar- 

Finally, field galaxies of deep surveys are fit- 
ted either by less massive "elliptical" scenarios of 
Mbar= 10^-^ *° " Mq or by "spiral" scenarios of 
Mbar=10^° *° Mq. The Hubble K diagram is alone un- 
able to separate the two models. In particular at low z 
(< 0.3), the predictions of the most massive elliptical and 
spiral models are similar. 

The K-z sequences are predicted with two other IMFs. 
The variation of results computed with either Rana & 
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Fig. 4. K-z sequences for Afbar = 10^ to 10^^ M© compared to field and radio galaxy samples in the Hubble iiT-band 
diagram. The IMF is from Rana & Basu, 1992, except for Afbar = 10^^ M©, the top-heavy IMF is also used (dotted 
line). The 10^^ Mq elliptical sequence (a) coincides with the K-z limit, up to z > 4. The spiral galaxy scenarios (b) are 
unable to fit distant radio galaxies at high z except for Mbar = 10^^ Mq, this value is however unacceptable at low z. 



Basu's IMF or the "top-heavy" (TH) IMF is presented 
on the Mbar = 10^^ M© sequence in Fig. g^b. Fig. El also 
shows the Mbar — M© sequence with the Salpeter 
IMF compared to Rana & Basu. The IMF effect on the K 
luminosity is less than 10% at all z. 

The Mbar, max scqueuce limit for massive ellipticals de- 
pends little on the formation redshift Zfor- Changing Zfor 
(>6.6, see Hu et al., 2002; Kodaira et al., 2003) in ellip- 
tical scenarios modifies the K-z relation only during the 
first Gyr (see Fig. 0Ji). 

The standard cosmology (Ho=65 km. s^^. Mpc^^, 
r2A/=0.3, riA—O.l) has been adopted. However, to study 
the sensitivity to the cosmological model, we compare the 
iiT-band observations to evolution models in the z (and not 
log z as previously shown) scale (Fig. |^ for various uni- 
verses. A flat universe without the cosmological constant 
(Ho=65 km. s~^. Mpc~^, ilM = l, ^Ia^O) is not compati- 
ble with the data. Both standard and open universes can 
fit the data; the variation of A/bar, max given by the two 
models is lower than 10 %. 

Fig. El is similar to Fig. 2|ib and Fig.|31but in the red- 
shift scale and includes more data: the catalogues 7C-I and 
7C-II (Wniott et al. 2003), 7C-III (Lacy et al. 2000) and 
6C* (Jarvis et al. 2001) are added to the De Breuck com- 
pilation. All our previous conclusions, in particular that 
the more powerful radio galaxy hosts are the more massive 
galaxies, up to 10^^ M©, are confirmed by these data. 

Emission lines from gas photoionization by massive 
stars are computed in our models. Fig. [7| shows pre- 
dictions for the main emission lines (Ha, [OIIIJA5007, 
[OIIJA3727), redshifted in the /f-band and compared to 
the data. The K-hand luminosity variation AK=Knoiine — 
Klines, due to ionization by massive stars, does not exceed 
0.25 mag in spirals, as in ellipticals. However, observa- 
tions show stronger emission lines than model predictions: 
other sources of ionization (shocks and photoionization by 



AGN) are known to be highly efficient in radio galaxies. 
From the powerful NIRC and the 7C radio sources (Fig. 
13), the AK/K relative excesses due to emission lines are 
<5% at high z. 

5. The K-z relation and the fragmentation limit 

The uniform galaxy mass limit Afbar,max= lO^^M© at all 
redshifts from the K-hand luminosities is striking. The 
critical mass of fragmentation for a self-gravitating cloud 
has been estimated by Rees & Ostriker (1977). The dy- 
namical model predicts that the kinetic energy of infall 
is first thermalized via shocks; further evolution then de- 
pends on the relative cooling tcooi and free-fall tff time- 
scales. High masses in the range 10^° to lO^^M© , under- 
going an efficient cooling, collapse at about the free-fall 
rate, fragment and possibly form stars while larger masses 
may experience a quasi-static contraction phase; they go 
into free fall only when they reach critical values of ra- 
dius and mass. The authors evaluated the critical mass 
Merit — lO^^M© between the two regimes that define 
the fragmentation limit. They also estimated the mass- 
independent critical size Vcrit —75 Kpc. 

The maximum baryonic mass Afbar,max delimitating 
the observed cut off of the galaxy distribution in the 
Hubble diagram is then similar to the theoretical fragmen- 
tation limit of galaxies. This clarifying explanation of the 
K-z relation is very simple. Moreover the corresponding 
critical size rcrit — 75 Kpc remarkably fits the typical di- 
ameter of Lya haloes (150 Kpc) surrounding active radio 
galaxies. The value Afbar,max at z ~0 is also dynamically 
confirmed by measurements of masses of nearby galaxies. 

This correspondance between the various mass esti- 
mates emphasizes the robustness of mass determination 
using K-hand stellar luminosities. This method of mass 
determination is thus potentially appropriate to all deep 
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Fig. 5. The A/bai— 10^'^ sequence shows the comparison of Rana & Basu IMF (long-dashed line) with Salpeter 
IMF(dotted line). The rest is similar to Fig.^. 



galaxy surveys which may improve our knowledge of mass 
evolution at high z. 

The Mbar,max Hmit for radio galaxies is also estimated 
from observations, up to z > 4 (Fig. implying strong 
constraints on galaxy evolution models. The most mas- 
sive galaxies which reach the fragmentation limit at z ~4 
formed their stellar mass within IGyr. This time-scale is 
too short to form the stellar mass by the hierarchical merg- 
ing of initial ~ IO^Mq building blocks as predicted by 
the CDM model. However these objects might be rare 
and their mode of formation exceptional. On the other 
hand, because distant radio galaxies are known to be in a 
dense medium of already evolved and luminous (masses 
~ 1O^^M0?) galaxy companions, the self-gravitational 
model may efficiently accumulate several lO^^M© within 
1 Gyr if the cooling process is highly efficient: the dissipa- 
tion by strong emission lines of the huge (diameter of ~ 
150 Kpc) ionized cocoon of the radio galaxies may explain 
the cooling efficiency at high redshifts. Chokshi (1997) al- 
ready evoked the Rees & Ostriker (1977) models to inter- 
pret distant radio galaxies, relating them to large density 
fluctuations undergoing isothermal free-fall and regulated 
by cooling H and He line emission. Her interpretation sug- 
gested intense stellar formation triggered by the radio jet. 
The present analysis of the Hubble ii'-diagram that con- 
firms large stellar populations in distant radio galaxies is 
unable to affirm that massive star formation was triggered 
by the radio jet, in the AGN environment or outwards 
from the overpressured cocoon. Only high spatial resolu- 
tion on extended radio galaxies will help to clarify such 
hypotheses. 

6. Masses of radio galaxy hosts 

In the literature, the stellar masses of radio galaxy hosts 
are tentatively related to their radio power: as an exam- 
ple, spheroids of stellar masses greater than 4 lO^^M© are 
more frequently observed in radio-loud than in radio-quiet 



galaxies (McLure & Dunlop, 2002). Willott et al. (2003) 
found systematic differences between the X-band magni- 
tude distributions in the Hubble diagram of the 7C, 6CE, 
6C* and SCR catalogues of increasing radio-power. The 
7C distribution is 0.55 mag fainter than the SCR at all red- 
shifts. Note that such a if-band magnitude difference at 
all redshifts means in our interpretation an average mass 
difference between the two catalogues. Fig|Slshows the re- 
spective mass distributions of the SCR catalogue (empty 
histogram) and of the 7C catalogue (grey line histogram) 
computed with our models. The two distributions are lim- 
ited by the mass limit Mbar,max — 10^^ Mq, confirming 
that the limit is independent of radio power. The more 
radio-powerful SCR catalogue is on average more massive 
(11.7 < Log (Mbar/MQ)<12.1) than the less powerful 7C 
catalogue (11.5 < Log (Mbar/Mo)< 12.0). The respective 
peaks of the two distributions differ by 0.2, corresponding 
to a variation of AMbar — 2 10^^ M0. The variation is not 
significantly sensitive to cosmological or IMF effects. 

The stellar and baryonic mass estimates from if-band 
luminosities with the help of evolution models is limited 
by uncertainties. According to Fig. the predicted se- 
quences computed with the open and standard universes 
are compatible with observations within 10% at the high- 
est z. As previously analyzed, the stellar masses of ellip- 
ticals vary within 10% depending on the adopted IMF. 
Another source of uncertainty at high redshifts is the cali- 
bration at z—Q. Depending on templates M\{Q,t{Q)) that 
are either the brightest cluster members or averaged el- 
liptical galaxies, the uncertainty may reach 10%. Strong 
emission lines at precise redshifts are seen in the Hubble 
K diagram so that they require a special fitting procedure. 
More uncertain is the nebular continuum of radio sources 
observed in the blue band (rest frame), not modelled in el- 
liptical galaxies. On the basis of starburst models, it may 
be lower than 20% at z >4. The global uncertainty on the 
upper limit of baryonic mass of galaxies given by mod- 
els is ~ 50%. That implies Mbar,max= 1 ± 0.5 IO^^Mq. 
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Fig. 6. Interpretation of the K-z diagram with models of various baryonic masses as in Fig. 01 or Fig but z is in 
a hnear scale. The catalogues 7C-I and 7C-II (Willott et al. 2003), 7C-III (Lacy et al. 2000) and 6C* (Jarvis et al. 
2001), are added to the sample from De Breuck et al., (2002). The sensitivity to cosmological parameters is analyzed 
on the 10^^ Mq sequence of elliptical galaxies (red full lines) for 3 cosmologies (open, flat and standard). Only open 
and standard (flat with a cosmological constant) universes are acceptable by data. 



For less luminous galaxies, the accuracy will depend on 
evolutionary scenarios by types. 



7. Discussion 

The main results obtained in this analysis are the follow- 
ing: first, the K-z relation is due to the galaxy mass limit 
lO-'^^M©; second, fitting the K-z sequence corresponds to 
the fragmentation limit predicted by gravitational mod- 
els; finally, we are able to identify massive galaxies of 
10^^ M0 at redshifts z ~4. The mass limit from Rees & 
Ostriker, (1977) based on dynamical arguments validates 
the baryonic mass determination based on the if-band 
luminosity. 

— Estimates are obtained using star formation scenar- 
ios with infall from the baryonic halo, winds, active 
and passive stellar evolutions. The star formation rate 
is proportional to the gas amount. We confirm that 
only the star formation scenarios of elliptical galax- 
ies are able to reproduce the stellar populations of 



the most powerful radio galaxy hosts. The robust- 
ness of our model PEGASE is to predict colors and 
SEDs of the typical observed templates at z = 0. 
Statistical values of such colors were estimated after 
corrections (inclination, extinction, aperture) by Fioc 
& Rocca-Volmerange (1999b). Scenarios of star forma- 
tion, galactic winds, infall, also respecting metal and 
dust evolution constraints in PEGASE, allow us to pre- 
dict the details of the evolution of the various mass 
components (FigUJ. 

Previous studies already tested the validity of the 
scenarios at higher z. The U to K multi-spectral 
faint galaxy counts from the HDF-N and Hawaii sur- 
veys were fitted with these scenarios (Fioc & Rocca- 
Volmerange 1999a). Another constraint of the valid- 
ity of the scenarios is given by the photometric red- 
shifts computed with the code Z-PEG. Used as high 
z templates, evolutionary scenarios give photometric 
redshifts close to spectroscopic redshifts with cr < 0.1 
(Le Borgne & Rocca-Volmerange 2002). 
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Fig. 7. Predictions through the if -band of the main emission hnes Ha, [OIII]A5007, [OII]A3727 from photoionization 
by massive stars for spirals and ellipticals. Comparison data are the SCR, 6C and 7C catalogues of radio sources 
(Willott et al. 2003). Galaxy ages depend on rcdshifts defined by fines (z=2.3, 3.3 and 5 for H^, [QUI] and [Oil]). The 
excess between observations and models is due to other ionization sources (shocks, AGN). 
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Fig. 8. Baryonic masses of 3C and 7C galaxies estimated from the elliptical scenario. The more radio-powerful 3CR 
galaxies (empty histogram) are on average more massive than the less radio-powerful 7C galaxies (grey line histogram). 
The difference D[Log(Mbar/ M©)] ~ 0.2 corresponds to AMbar =i 2 IO^^Mq. 



Does the observed dispersion of the K-z relation 
weaken our results? One fundamental property of the 
relation is its uniformity (Lilly & Longair, 1984) with 
an increase of the dispersion at redshifts z > 2 (Ealcs 
et al, 1997; Lacy et al., 2000). McCarthy et al. (1999) 
compared the present samples (SCR, 6C and NIRC) 
to the Molonglo (MRC/lJy) sample, complete in the 
sense of containing > 99% of all sources above its flux 
limit (S408MHZ > 0.95 Jy). The authors conclude to 
the consistency of the SCR, MRC and 6C samples 
within 1(7 at z < 0.5. At higher z, Willott et al.(200S) 
combine the 6C, 6C* samples and the 7C sample, com- 



plete above the flux limit (Si5iMffz > 0.5 Jy) in three 
regions of the sky. They confirm the homogeneity of 
the stellar populations formed at high z and follow- 
ing a passive evolution. The variation from z=0.05 to 
3 is lower than 5%, in agreement with Jarvis et al. 
(2001). They also found that the interpretation by us- 
ing two models (no-cvohition and zjor=5 ) are not in 
agreement with the data. We do not study them in the 
present paper. Their conclusion is that the decrease of 
radio galaxy masses at high redshifts is small, presum- 
ing a very rapid evolution for radio sources at high 
z. 
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Regarding the debate on the IMF as the explanation 
for divergences between observables and models, we 
confirm that K magnitudes are faintly sensitive to the 
IMF (Tinsley, 1972). However a reasonable value of 
the lower mass is essential since the bulk of stars ap- 
pearing in the iC-band at z—Q is made up of low mass 
stars. The value of 0.1 Mq, adopted for all IMFs, is 
compatible with the data and with theoretical star 
formation principles. By respecting this constraint on 
the lower star limit, changing the IMF (including top- 
heavy) slope will not significantly modify our conclu- 
sions within the error bars. 

Remarks about the comparison of dissipative collapse 
and hierarchical merging scenarios are still prelimi- 
nary. However, some general ideas may be derived from 
our conclusions. Stellar masses predicted aX z = for 
ellipticals and spirals (respectively 55 % and > 45 % of 
the initial cloud of mass Mbar) are in agreement with 
typical values for nearby field galaxies. They range 
from 10^ M0 for dwarf to 10^^ M© for giant galax- 
ies. Field galaxies correspond to long time-scales (> 10 
Gyr) of less massive clouds, which may reach by hier- 
archical merging a final mass of 10^^ M©, the fragmen- 
tation limit. The time-scales of mass accumulation, gas 
depletion and star formation of massive ellipticals at 
high z are, however, incompatible with building blocks 
of 10^ Mq, proposed by the CDM model. Our results 
cannot exclude some merging process of massive (a 
few 10^^ Mq) blocks in the dense environment of radio 
galaxies or proto-clusters. 

Is the presence of black holes in radio sources linked to 
host elliptical formation? To interpret massive ellipti- 
cal galaxies at high z requires the accumulation of the 
dynamical mass on a short time-scale. It was suggested 
that the embedded black hole is associated with a rapid 
mass accumulation and star formation process. The 
more massive the black hole, the deeper the potential 
wall; this does not necessarily favor short time-scales 
for accumulation of baryonic and stellar mass in the 
environment of the potential wall. In the absence of a 
magnetic flux, a massive cloud evolves as a single unit, 
after having shed its angular momentum via instabili- 
ties. The dense environment may favor the apparition 
of a black hole (Silk & Rees, 1998). The high density of 
galaxy companions within the ionized halos of distant 
radio galaxies (of critical size 150 kpc) has been shown 
by integral field spectroscopy (Rocca-Volmerange et 
al., 1994), and narrow-band imaging of the Lya line 
(Kurk et al., 2003). Moreover the interactive role of 
a radio jet to trigger star formation is still a debated 
subject. Such an hypothesis is favored by the observa- 
tional relations between black hole masses and abso- 
lute blue luminosities or velocity dispersion (Ferrarese 
& Merritt, 2000; Gebhardt et al. 2000; Magorrian et 
al. 1988). However it is too early to identify the nature 
of the relation between stellar and black hole masses 
(McLure & Dunlop, 2002). The bulk of star formation 
could happen during the free-fall time while the ac- 



tive nucleus is fueled. The black hole and the stellar 
mass would then simultaneously grow; details of such 
processes would, however, need clarification. 
Recent deep X-ray surveys which follow the evolution 
of accretion onto supermassive black holes confirm two 
modes of accretion and of black hole growth. While the 
rare, high luminosity QSOs could form rather early in 
the universe, a late evolution of low-luminosity Seyfert 
populations is required by observations (Hasinger et 
al., 2003). Because high luminosity QSOs are embed- 
ded in massive elliptical galaxies while low-luminosity 
Seyfert galaxies are preferentially found in spiral galax- 
ies, the two modes of star formation derived from the 
i^T-band Hubble diagram can be related to the two ac- 
cretion modes derived from X-rays. Moreover the K- 
band Hubble diagram of sub-mm and hyperluminous 
galaxies (Serjeant et al. , 2003) recently confirms that 
the location of galaxies on the K-z relation may be 
related to the presence of the most massive AGNs, in 
perfect agreement with our conclusions. 

A new factor of uncertainty could be the ultra-violet 
continuum of radio sources. However because galaxies are 
observed in the K band, only radio galaxies at z > 4 will 
be affected by this effect. 

Complementary constraints are awaited from im- 
proved spatial resolution (morphology, surface brightness, 
velocity fields). The (1 -I- z)^ surface brightness fading bi- 
ases high z observations towards bulges rather than disks. 
This may explain why, at z ~ 4, field galaxies appear 
on lower Mbar sequences (10^ to 10^^-^Mq) than radio 
galaxies of high surface brightnesses. We may also evoke 
the high density of galaxies within their diameter-limited 
sample to explain massive IO^^Mq radio sources at high 
redshift. Only a better spatial resolution, as given by in- 
tegral field units on the new generation of telescopes, will 
clarify the spatial confusion of the central radio galaxy 
and its galaxy companions. 

The baryonic mass predicted by our modeling of lumi- 
nosity evolution is a minimal value. From Fig. |21at z ~0 
(14 Gyr), a mass fraction 1—Mgai (35 % for ellipticals, 10% 
for spirals) is still in the halo. The galaxy halo could also 
contain another baryonic component (such as hot gas), 
not taken into account in our mass budget. 

8. Conclusion 

The main results of this analysis are the following: 

i) Stellar masses of galaxies, and the corresponding 
baryonic masses of their progenitor clouds, are estimated 
from the Lk luminosities at any z using robust scenarios 
of star formation evolution. 

ii) The brightest luminosity limit of the Hubble K di- 
agram corresponds to the most massive elliptical galaxies 
of baryonic masses Mbar, max — 10^^ Mq (about 55 % is 
in stellar mass at z—Q). These are the typical hosts of 
powerful radio galaxies. 
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iii) The maximum limit of baryonic mass 
Mbar,max — 10^^ M0 , estimated from the stellar lu- 
minosity in the ii'-band, corresponds to the dynamical 
parameter M^rit, the fragmentation hmit of the cooling- 
gravitation balance of a self-gravitating cloud. This 
confirms the validity of the mass estimates from evolution 
scenarios. 

iv) Radio galaxies of 10^^ Mq are found at z > 4 with 
star formation and mass accumulation time-scales shorter 
than IGyr. The radius of the ionized cocoon is about 
'''crit =75kpc. Such extreme time-scales in powerful radio 
galaxies suggest that the formation of a massive black hole 
embedded in the host galaxy is not independent of the star 
formation process. 

v) The choice of the IMF slope is not crucial but the 
lower star mass needs to be 0.1 to form the bulk of low 
mass stars observed in templates of elliptical galaxies at 
z=0. 

vi) Galaxy deep surveys are covered by sequences of 
lower baryonic masses from Mbar= 10^ to 10^^ in 
which the evolution by hierarchical merging, if it is proven, 
remains however mass limited by fragmentation. 

For the future, we will need to apply such a method 
of mass determination to the galaxy populations of large 
surveys, in particular discovered either by SIRTF or by 
the ground-based 10m telescopes in the stellar infra-red 
to follow the evolution processes of galaxies. 
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